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ABSTRACT 

We report on the detection in Sloan Digital Sky Survey data of at least three, 
roughly parallel components in a 65°-long, stellar stream complex previously 
identified with the Anticenter or Monoceros Ring. The three-stream complex 
varies in width from 4° to 6° along its length and appears to be made up of 
two or more narrow substreams as well as a broader, diffuse component. The 
width and complexity of the stream indicate that the progenitor was likely a 
dwarf galaxy of significant size and mass. The stream is 8.9 kpc distant and is 
oriented almost perpendicularly to our line of sight. The visible portion of the 
stream does not pass near any known dwarf galaxies and a preliminary orbit 
does not point to any viable progenitor candidates. Orbits for the narrower 
substreams can be modeled with velocity offsets from the broad component of 
8 km s~^. We suggest that the broad component is likely to be the remains of a 
dwarf galaxy, while the narrower streams constitute the remnants of dynamically 
distinct components which may have included a native population of globular 
clusters. While the color of the main sequence turn-off is not unlike that for 
the Monoceros Ring, neither the visible stream nor any reasonable projection of 
its orbit passes through Monoceros or Canis Major, and we conclude that this 
stream is probably unrelated to the overdensities found in these regions. 

Subject headings: Galaxy: Structure — Galaxy: Halo 



1. Introduction 

The value of large scale digital sky surveys to studies of Galactic structure and the Local 
Group has become abundantly clear in recent years, and particularly in the last few months. 
In addition to the large scale features attributed to past galaxy accretion events (Yanny et 
al. 2003; Majewski et al. 2003; Rocha-Pinto et al. 2004), Sloan Digital Sky Survey (SDSS) 
data were used to detect the remarkably strong tidal tails of Palomar 5 (Odenkirchen et al. 



- 2 - 



2001; Rockosi et al. 2002; Odenkirchen et al. 2003; Grillmair & Dionatos 2006a), NGC 
5466 (Belokurov et al. 2006a; Grillmair & Johnson 2006), and streams due to extant or 
extinct globular clusters (Grillmair & Dionatos 2006b) and dwarf galaxies (Belokurov et al. 
2006b; Grillmair 2006). Willman et al. (2005), Zucker et al. (2006), and Belokurov et al. 
(2006c) recently used SDSS data to discover several new dwarf satellites of the Milky Way. 

In this paper we continue our analysis of the SDSS database to search for extended 
structures in the Galactic halo. We briefly describe our analysis in Section 2. We discuss 
the structure of a stream complex, attributed to the Monoceros Ring by Belokurov et al. 
(2006b) in Section 3, estimate distances in Section 3.1, and put constraints on the orbit in 
Section 3.2. 



2. Data Analysis 

Data comprising ugriz photometry for stars in the region 108° < a < 162° and —4° < 
6 < 68° were extracted from Data Release 5 of the SDSS using the SDSS CasJobs query 
system. The data were analyzed using the matched filter technique employed by Grillmair 
& Johnson (2006), Grillmair & Dionatos (2006a), Grillmair & Dionatos (2006b), Grillmair 
(2006) and described in detail by Rockosi et al. (2002). This technique is made necessary 
by the fact that, over the magnitude range and over the region of sky we are considering, the 
foreground disk stars outnumber the more distant stars in the Galactic halo by some three 
orders of magnitude. Applied in the color-magnitude domain, the matched filter is a means 
by which we can optimally differentiate between two populations. 

We used the SDSS photometry to create a color-magnitude density or Hess diagram for 
both stars of interest and for the foreground population. Dividing the former by the latter, 
we generated an array of relative weights which we used as an optimal color-magnitude filter. 
We generated the search filter using the color-magnitude distribution of stars in M 13. Owing 
to M 13's proximity, we have a better measure of its luminosity function and the effects of 
SDSS completeness than we do for other objects in the DR5 field. We used all stars with 
lb < g < 22.5. We dereddended the SDSS photometry as a function of position on the sky 
using the DIRBE/IRAS dust maps of Schlegel, Finkbeiner, & Davis (1998). A single Hess 
diagram for field stars was generated using 1.2 x 10^ stars spread over ~ 2200 deg^ of DR5. 
We applied the M 13 filter to the entire survey area, and the resulting weighted star counts 
were summed by location on the sky to produce a two dimensional probability map. 

In Figure 1 we show the final filtered star count distribution, using a filter matched to 
the g,g — i color magnitude distribution of stars in M 13 but shifted faintwards by 0.3 mags. 
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An exponential surface function has been subtracted from the image to remove the rapid 
rise in the number of disk stars at low Galactic latitudes, and a 4th-order polynomial surface 
has been used to remove remaining large scale gradients. The final image has been binned 
to a pixel size of 0.08° and smoothed using a Gaussian kernel with a — 0.2°. 

3. Discussion 

Apparent in Figure 1 is a broad and complex stream running from north to south across 
the field. Portions of this stream are visible in Ncwberg et al. (2002) and Bclokurov et al. 
(2006b). The narrow, curved stream running towards the northeast is part of the 63°-long 
globular cluster stream found by Grillmair & Dionatos (2006b), hereafter referred to as 
GD-1. The broad, cast-west stream just above the main gap in the data is the Sagittarius 
stream discussed by Belokurov et al. (2006b), though somewhat muted by our filtering 
due to its greater distance. The stream complex of interest extends from (R.A., decl.) = 
(126.4°, -0.7°) to (R.A., decl.) = (133.9°, 64.2°), and runs in a 65°, nearly great circle path 
from Ursa Major in the north to Hydra in the south. The stream is truncated at both the 
southern and northern ends by the limits of the available data. We note that there is also an 
apparent concentration of stars at (R.A., decl.) = (134°, 3.4°), surrounded by faint, banded 
substructure roughly parallel to that of the main western stream. 

Prom the reddening map of Schlegel, Pinkbeiner, & Davis (1998), the maximum values 
of E{B — y) are 0.2 (near the northern tip), with typical values near 0.03 along the 
remainder of the stream. There arc diminutions here and there in the stream that could be 
attributed to regions of higher reddening, but there are no long features with a north-south 
orientation which could be held to account for either the appearance of the stream as a 
whole, or for the different components within it. 

Sampling at several representative points, the main stream complex appears to be about 
5° wide on average. This is significantly broader than the globular cluster streams found by 
Odcnkirchcn et al. (2003); Grillmair & Johnson (2006); Grillmair & Dionatos (2006a) and 
Grillmair & Dionatos (2006b). 5° corresponds to about 800 pc at our estimated distance to 
the stream (see below), which is much larger than the tidal diameters of globular clusters. 
We conclude that the progenitor was considerably more extended than a globular cluster 
and was most likely a dwarf galaxy. 

The stream is clearly not just a broad swath of stars. Rather, it appears to be made up 
of a ~ 2° wide broad component running roughly along the center of the stream complex, 
and at least two narrower, l°-wide streams to the east and west of the broad component. 
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The fine structure in the stream complex is illustrated in Figure 2, where we have made 
several east-west slices across the complex at various declinations. Comparing Figures 1 and 
2, there are indications of additional structure within the broad component, and still other, 
more tenuous parallel streams to the east and west of the three major components. 

Integrating the background subtracted, weighted star counts along the stream over a 
width of ~ 5° wc find that the total number of stars in the discernible stream down to 
g = 22.5 is 9200 ± 1500. The mean surface density of stars in the southern portion of the 
stream ( — 1° < S < 14°) is about 67 stars deg^^, which is roughly twice the 32 stars deg^^ 
(corrected for cos (5)) found in the northern section (39° < S < 52°). The highest local 
surface densities exceed 200 stars deg^'^. If we assume a globular cluster- like luminosity 
function in the stream, then we can use the color transformation equations of Smith et al. 
(2002), the deep M 4 luminosity function of Richer et al. (2002), and the mass-luminosity 
relation of ? to extrapolate to fainter magnitudes. Integrating over 3 < My < 17, we 
estimate a total number of stars in the visible portion of the stream of 33, 000 ±5400, a total 
luminosity of 5500 ± DOOL©, and a total estimated mass of 9300 ± 15OOM0. 

3.1. Color-Magnitude Distribution and Distance to the Stream 

In Figure 3 we show g,g — i color- magnitude distributions for the stream stars, extracted 
by generating a Hess diagrams of stars lying along 2°-wide regions covering the western 
and eastern halves of the stream complex and subtracting a similar field star distribution 
sampled over ^ 500 deg^ to the east and west of the stream complex. Despite the somewhat 
limited statistics, a clear signature of the turn-off and main sequence is evident in the stream 
population. Moreover, the turn-off regions of the distributions match the dcrcddcned, shifted 
main sequence locus of M 13 fairly well. The turn-offs in the eastern and western halves of 
the complex lie at dereddened g — i — 0.27 and 0.3, respectively. A similar diagram in g — r 
yields dereddened turn-off colors oi g — r — 0.23 and 0.25, respectively, with the difference 
most likely a due to ~ 0.02 mag measurement uncertainty. These estimates lie between 
values of 0.26 and 0.22 measured by Newberg et al. (2002) for the Monoceros Ring and the 
Sagittarius stream, respectively. They are also much bluer than Newberg et al. (2002) 's 
estimate of dereddened g — r = 0.40 for thick disk stars. Based solely on turn-off color, 
we can therefore rule out association of the stream complex with thick disk stars. However, 
within the uncertainties, the turn-off colors are consistent with those of either the Monoceros 
Ring or the Sagittarius stream. 

Varying the magnitude shift applied to M 13's main sequence locus from -1.0 to +3.0 
mag, we measured the foreground-subtracted, mean surface density of stream stars in the 
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regions -1° < 5 < 9°, 17° < 5 < 39°, and 39° < 5 < 63°. To avoid potential problems 
related to a difference in age between M 13 and the stream stars, we used only the portion 
of the filter with 19.5 < g < 22.5, where the bright cutoff is 0.8 mag below M 13's main 
sequence turn-off. Though this reduces the stream contrast somewhat, it provides sufficient 
integrated signal-to-noise to enable main sequence fitting. 

Fitting Gaussians to the mean surface densities as a function of magnitude shift (e.g. 
Grillmair & Dionatos (2006b), we find that the highest contrasts occur for shifts of +0.31, 
+0.29, and +0.37 mag for the eastern half of the stream complex over the declination ranges 
given above, respectively. For the western half, we find that the filter response peaks at 
shifts of +0.34, +0.25, and +0.38 mag, respectively. The magnitude shifts for the two halves 
are highly consistent with one another, and we conclude that there is no significant distance 
offset from one side of the complex to the other. Adopting a distance to M 13 of 7.7 kpc 
(Harris 1996) we find an average hehocentric distance oi d — 8.9 ± 0.2 kpc. The stream is 
roughly perpendicular to our line of sight and slightly curved about the Galactic center as 
expected. Our distance estimate is in excellent agreement with the ^ 9 kpc found by Ibata 
et al. (2003) in their WFS-0801 field, which is situated on the western edge of the stream 
complex at a = 120.5, 6 = 40.3. 



3.2. Constraints on the Orbit 

The visible portion of the stream complex spans the Galactic anticenter direction and, 
projecting a great circle path, is inclined by 35° to the Galactic plane. Though we are 
currently limited by a lack of velocity information, for a given model of the Galactic poten- 
tial the progenitor's orbit is actually fairly well constrained by the observed distance and 
orientation of the stream. Using the Galactic model of Allen & Santillan (1991) (which 
includes a disk, bulge, and spherical halo, and which Grillmair & Johnson (2006); Grillmair 
& Dionatos (2006a) and Grillmair & Dionatos (2006b) found to work reasonably well for 
NGC 5466, Pal 5, and GD-1), we use a least squares method to fit both the orientation on 
the sky and the distance measurements in Section 3.1. In addition to a number of normal 
points lying along the central component of the stream, we chose as a velocity fiducial point 
a position at the northern end of the stream at (R.A., decl) = (125.463°, 51.492°). 

If we allow the proper motions to be free ranging and uninteresting parameters, the 
model which best fits the data predicts vlsr = — 18 ± 10 km s~^ at the fiducial point, where 
the uncertainty corresponds to the 95% confidence interval. A projection of this orbit is 
labeled C in Figure 4. We note that the uncertainty is primarily determined by the large 
lever arm over which it has been possible to measure relative distances. The 95% range in 
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vlsr in turn predicts a range in perigalactic and apogalactic radii of 6.6 < Rp < 6.9 kpc and 
16.8 < Ra < 17.3 kpc. Of course, these ranges do not take into account uncertainties in the 
absolute distance of the stream (which depends on the uncertainty in M 13's distance) or of 
the vahdity of Allen & Santillan (1991) 's Galactic model. 

Given the very similar distances estimated for the eastern and western portions of the 
stream, it is highly unlikely that the complex could be a superposition of multiple wraps 
around the Galaxy. The two narrower streams (E and W in Figure 4) can be reasonably 
well modeled by 0.18 mas yr^^ offsets in east-west proper motion. At the distance of the 
stream complex this amounts to ~ 8 km s~^. The entire stream complex is therefore likely 
to be the remains of a dwarf galaxy which contained distinct components spanning a range 
of binding energies. The narrower streams might, for example, be the remnants of the parent 
galaxy's globular cluster population. Piecing together the original structure and evolution 
of the stream's progenitor will require detailed N-body modeling. 

Integrating orbits for parameter sets spanning the range above, we find that, with the 
exception of the Sagittarius dE, there are no known dwarf galaxies within 5° of the projected 
orbit. The Sagittarius dE hes 3.1° from the projection of the best-fit orbit, but the orbital 
planes of the the Sagittarius dE and the new stream are clearly distinct (Figure 4); we 
attribute this apparent proximity to the expected confluence of orbit projections in the 
direction of the Galactic center and not to any physical association between them. 

The orientation of the stream on the sky puts fairly strict limits on the plane of the 
orbit. The visible portion of the stream passes through Lynx, Cancer, and Hydra. We 
find no reasonable combination of parameters that would place the southern projection of 
the stream in either Monoceros or Canis Major. Nor docs the visible extent of the stream 
complex fit either the prograde or retrograde models of the Monoceros stream computed by 
Penarrubia et al. (2005). Thus, even while the turn-off colors appear to be similar, we 
conclude on orbital grounds that this stream complex is unlikely to be related to cither the 
Monoceros stream or the Canis Major overdensity. 

Funding for the creation and distribution of the SDSS Archive has been provided by 
the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics and 
Space Administration, the National Science Foundation, the U.S. Department of Energy, the 
Japanese Monbukagakusho, and the Max Planck Society. 
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- 7- 



REFERENCES 

Allen, C, & Santillan, A. 1991, Rev. Mex. Astron. Astrofis., 22, 255 

Belokurov, V., Evans, N. W., Irwin, M. J., Hewett, P. C, & Wilkinson, M. I. 2005, ApJ, 
637, L29 

Belokurov, V. et al. 2006, ApJ, 642, 137 
Belokurov, V. et al. 2006, ApJ, 647, Llll 
Grillmair, C. J. 2006, ApJ, 645, 37 
Grillmair, C. J., & Johnson, R. 2006, ApJ, 639, L17 
Grillmair, C. J. & Dionatos, O. 2006b, ApJ, 641, L37 
Grillmair, C. J., & Dionatos, O. 2006, ApJ, 643, 17 
Harris, W. E. 1996, AJ, 112, 1487 

Ibata, R. A., Irwin, M. J., Lewis, G. F., Ferguson, A. M. N., & Tanvir, N. 2003, MNRAS, 
340, 21 

Majewski, S. R., Skrutskie, M. F., Weinberg, M. D., & Ostheimer, J. C. 2003, ApJ, 599, 
1082 

Newberg, H. J. et al. 2002, ApJ, 569, 245 

Odenkirchen, M., Grebel, E. K., Rockosi, C. M., Dehnen, W., Ibata, R., Rix, H., Solte, A., 
Wolf, C, Anderson, J. E. Jr., Bahcall, N. A, Brinkmann, J., Csabai, I., Hennessy, G., 
Hindsley, R. B., Ivezic, Z., Lupton, R., Munn, J. A., Pier, J. R., Stoughton, C, York, 
D. G. 2001, ApJ, 548, 1650 

Odenkirchen, M., Grebel, E. K., Dehnen, W., Rix, H., Yanny, B., Newberg, H., Rockosi, C. 
M., Martinez-Delgado, D., Brinkmann, J., Pier, J. R. 2003, AJ, 126, 2385 

Penarrubia, J., Martinez-Delgado, D., Rix, H. W., Gomez-Flechoso, M. A., Munn, J., New- 
berg, H., Bell, E. F., Yanny, B., Zucker, D., & Grebel, E. K. 2005, ApJ, 626, 128 

Richer, H. B. et al. 2002, ApJ, 574, L151 

Rocha-Pinto, H. J., Majewski, S. R., Skrutskie, M. F., Crane, J. D., Patterson, R. J. 2004, 
ApJ, 615, 732 



-8- 

Rockosi, C. M., Odenkirchen, M., Grebel, E. K., Dehnen, W., Cudworth, K. M., Gunn, J. 
E., York, D. G., Brinkmann, J., Hennessy, G. S., & Ivezic, Z. 2002, AJ, 124, 349 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Smith, J. A. et al. 2002, AJ, 123, 2121 

Willman, B. et al. 2005, AJ, 129, 2692 

Yanny, B. et al. 2003, ApJ, 588, 824 

Zucker, D. B. 2006, ApJ, submitted. 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-9- 




Fig. 1. — Smoothed, summed weight image of the SDSS field after subtraction of both an 
exponential and a 4th-order polynomial surface fit. Darker areas indicate higher surface 
densities. The weight image has been smoothed with a Gaussian kernel with a — 0.2°. The 
solid gray areas arc either data missing from DR5, or clusters or bright stars which have been 
masked out prior to analysis. The stream complex runs from (R.A., decl) = (126.4°, —0.7°) 
to approximately (R.A., decl.) = (133.9°, 64.2°). The feature at (R.A., decl.) = (143°, 30°) 
is a portion of the cold stellar stream discovered by Grillmair & Dionatos (2006b). 



-10- 




140 135 130 125 120 

R.A. (deg. J2000) 



Fig. 2. — Profiles across tlic stream complex at constant declination, after smoothing the 
weighted star counts with a Gaussian kernel with a — 0.25° and summing over 2.4°-wide 
bands in dechnation. E, C, and W indicate the Eastern, Central, and Western components. 
EBS denotes the Easter Banded Structure alluded to in the text. The profiles have been 
offset vertically for clarity. Note indications of still other streams, both within the central 
stream and to the east and west of the complex. 
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Fig. 3. — The dereddened, background subtracted, color-magnitude distribution of stream 
stars. The left-hand panel shows the C-M distribution of stars in the western portion of 
the stream complex, while the right-hand panel shows the distribution of stars in a 2°-wide 
strip on the eastern side of the complex. The main sequence and the turn-off are clearly 
distinguishable. The solid line shows the dereddened locus of giant branch and main sequence 
stars as measured in DR5 for M 13, shifted faintwards by 0.3 mags. 
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Fig. 4. — Contour plot of Figure 1, showing projections of the best fitting orbits to each of 
the eastern, central, and western stream components. GD-1 is the globular cluster stream 
found by Grillmair & Dionatos (2006b), and other designations are as in Figure 2. 



